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ABSTRACT 

IRAS 13224-3809 was observed in 2011 for 500 ks with the XMM-Newton observatory. 
We detect highly significant X-ray lags between soft (0.3 - 1 keV) and hard (1.2 - 5 
keV) energies. The hard band lags the soft at low frequencies (i.e. hard lag), while 
the opposite (i.e. soft lag) is observed at high frequencies. In this paper, we study 
the lag during flaring and quiescent periods. We find that the frequency and absolute 
amplitude of the soft lag is different during high-flux and low-flux periods. During 
the low flux intervals, the soft lag is detected at higher frequencies and with smaller 
amplitude. Assuming that the soft lag is associated with the light travel time between 
primary and reprocessed emission, this behaviour suggests that the X-ray source is 
more compact during low-flux intervals, and irradiates smaller radii of the accretion 
disc (likely because of light bending effects). We continue with an investigation of the 
lag dependence on energy, and find that isolating the low-flux periods reveals a strong 
lag signature at the Fe Ka line energy, similar to results found using 1.3 Ms of data 
on another well known Narrow-Line Seyfcrt I galaxy, 1H0707-495. 
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1 INTRODUCTION 

The X-ray emitting region in accreting black hole systems is 
not well understood. Generally, it is agreed that the angular 
momentum of accreting gas onto a black hole forms an op- 
tically thick, geometrically thin accreti on disc that radiates 
as a series of blackbody components (jShakura fc SunvaevI 
1 19731 ). For supermassive black holes, this thermal emis- 
sion peaks in the UV. Some of the thermal disc photons 
are then inverse Compton upscattered to X-ray energies by 
mildly relativistic electrons in a hot cloud, called the corona 
l|Haardt fc MaraschillT99ll . Il993l '). While the X-ray emission 
mechanism is relatively well understood, the geometry of the 
corona where the X-rays are produced is an area of active 
research. In this paper, we probe the geometry of the corona 
through an analysis of the X-ray reverberation lags. 

X-ray variability studies have been influential in start- 
ing to understand the size of the emitting region. AGN are 
known to be extremely variable in the X-ray band, com- 
monly doubling in amplitude on timescales of less than one 
day for radio-quiet sources (jMcHardvl 1 19881 ). Substantial 
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variability occurs on timescales of the light crossing time of 
the source (or greater) , and therefore we can calculate an ug- 
per li mit on the size of the emitting region, R < c5t ( Tcrrel]| 
Il967l ). From observations of fast variability in AGN, we know 
that the corona must be compact, likely within ~ 100 Vg. 
Recent observations from the X-ray variability seen from 
microlensed quasars al so show that the corona is compact 
to within 10 or l ess (|Chartas et al.ll2009l : Ibai et al.ll201Gl : 
IChartas et "al]|2012l ). 

Variability studies of black hole binaries led to the 
discovery of the X-ray lag, whe re usually hard pho- 
tons lag behind soft photons fe.g. iMivamoto et al.l Il988l : 
iNowak fc Vaughanll996l : 'Nowak et al."l999') . This effect was 
also later observed in AGN (Papadakis, Nandra, & Kazanai 
I2OOII : iMcHardv et al.l [20041 ). The origin of the hard lag is 
still not well understood, though in the prevailing model, 
the multiplicative propagation effects of fluctuations in the 
accretion flow modulate regio ns emitting soft phot ons before 
those emitting hard photons (|Kotov. Churazov.~fc Gilfanovl 
I2OOII : lArevalo fc Uttlevll2006l ). 

Recent developments in X-ray variability have revealed 
a new type of lag, where soft photons lag behind hard pho- 
tons. This lag, interpreted as a reverberation lag, offers a 
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new perspective in which to study the X-ray emitting re- 
gion in AGN. Understanding the phase lag between light 
curves of different energy bands gives us an orthogonal ap- 
proach to testing physical models that are often degenerate 
in the time-integrated energy spectrum. 

Reverberation lags were first discovered in 
Narrow-Line Seyfer t 1 (NLSl) galaxy, 1H0707-495 
l|Fabian et al.l |2009| ). Since then, lags have been 
detect ed more than a dozen other Seyfert galax- 
ies (|Emmanoulopoulos. McHardv. fc Papadakid I2OIII : 
iDe Marco et al.ll201ll . I2OI2I : IZoghbi et al.ll2012l). ir icludiii g 
most recently IRAS 13224-3809 (|Fabian etall l2012bl ). 
Reverberation lags have also been observed on the order 
of milliseconds or less in black hole binary GX 339-4 
l|Uttlev et al.ll201ll ). 

The soft lag s can be understood with the standard re- 
flection scenario l|Guilbert fc Reedl 19881 : iLightman fc White! 
Il988l ). In this model, the bulk of the X-ray continuum is pro- 
duced in the corona. Some continuum photons fall towards 
the surface of the accretion disc, are absorbed by photo- 
electric absorption, and reprocessed into an X-ray reflection 
spectrum c omprised of a reflecti on continuum and emission 
lines (e.g. IRoss fc Fabianllioosi ). When reflection arises in 
the accretion disc, the whole reflection spectrum is blurred 
due to Doppler and other relativistic effects close to the 
black hole. Soft lags, where the direct continuum variations 
lead those in the soft-band reflection, are then interpreted as 
the additional light travel time taken by the reflected paths 
between the corona and the inner accretion disc. 

IRAS 13224-3809 {z = 0.066) i s one of the most X-ray 
varia ble Seyfert 1 galaxies known (|Boller. Brandt, fc FinkI 
Il996l ). and therefore is a useful source with which to probe 
the environments of the innermost regions. It was first ob- 
serve d with XMM-New t on in 2002 for 64 ks ijBolIer et all 
I2OO3I : iGallo et al.l l2004l : iPonti et al.l I2OI0I I, and most re- 
cently, for 5 00 ks, which led to a significant detection of 
the soft lag (|Fabian et al.|[2012bh . In this letter, we follow 
up that recent work with an analysis of the soft lags of 
IRAS 13224-3809 in the low and high flux states. 



2 OBSERVATIONS AND DATA REDUCTION 



The XMM-Newton satellite (|jansen et al.l 120011 ) observed 
IRAS 13224-3809 for 500 ks over four orbits from 2011 
July 19 to 2011 July 29 (Obs. IDs 0673580101, 0673580201, 
0673580301, 0673580401). For this analysis of the reverbera- 
tion lags, we focus on the high time-re solution data from the 
EPIC-pn camera (IStriider et al.ll200ll ). The first observation 
was taken in full window imaging model, and next three in 
large window imaging mode. All of the data were reduced 
in the same way, using the XMM-Newton Science Analysis 
System (SAS v. 11. 0.0) and the newest calibrat ion files. The 
details of the data reduction are explained in iFabian et al.l 
l|2012b[ ). 

The data were cleaned for high background flares, which 
resulted in a total exposure time of ~ 300 ks. The data were 
selected with the condition pattern < 4. Pile-up effects 
were not significant in any of the observations. 

The source light curve was extracted from circular re- 
gions of radius 35 arcsec, which were centered on the max- 
imum source emission. The background light curves were 
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Figure 1. Lag-frequency spectrum for the 500 ks observation. 
The lag is calculated between the soft energy band (0.3 - 1. keV) 
and the hard band (1.2 - 5. keV). The most negative lag (at 
v = A.lx 10"* Hz) is -92 ± 31 s. 



chosen from a circular region of the same size, and were 
the same distance to the readout node as the source region. 
The background subtracted light curves were produced us- 
ing the tool EPICLCCORR. The light curves ranged in length 
from 83 ks to 124 ks with 10 s bins. 



3 RESULTS 

3.1 Lag vs. Frequency 

The lag is co mputed in the usual way following the formulae 
described in lNowak et al.l (|l999l ). The frequency-dependent 
lag is calculated between light curves in the soft (0.3 - 
1. keV) and hard (1.2 - 5 keV) energy bands, which are 
dominated by the soft excess (likely disc-refiec tion) and the 
primary powerlaw-fike e mission, respectively (iFabian et al.l 
l2012bl : [Ponti et al.llioTol '). The Fourier transform of a light 
curve contains an amplitude and complex phase, such that 
the Fourier transform of soft light curve s{t) is written 
S = ISIe^^^. The cross product of the soft and hard band 
Fourier transforms, H* S = |i?||5|e'('*»"'^''\ provides the 
phase difference between the two energy bands. This phase 
difference is converted back into a frequency-dependent time 
lag, where r(/) = (</)s — <j)h)/'2nf- Given this definition for 
the lag, a negative amplitude lag means that the soft band 
lags behind the hard band, while the positive amplitude lag 
shows the hard band lagging. We take note that the phase 
lag is defined over the interva l (— tt.tt), which causes p hase 
wrapping at high frequencies (|Nowak fc Vaughanlll99q ). 

Fig. [T]shows frequenc y-dependent t i me lag between the 
soft and hard bands, from IFabian et aU (|2012bl ). At low fre- 
quencies, we find a large positive lag where the hard fiux lags 
behind the soft ffux by hundreds of seconds. At higher fre- 
quencies, u ~ [2-10] X 10"" Hz, there is a negative lag where 
the soft ffux lags the hard by ~ 100 s. Above ~ 1.5 x 10"^ 
Hz, the signal becomes dominated by Poisson noise. 



3.2 Flux-resolved analysis 

IRAS 13224-3809 is a highly variable source that exhibited 
distinct ffaring and quiescent periods during the 500 ks ob- 
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Figure 2. Hard band (1.2 — 5 keV) and soft band (0.3 - 1 keV) light curves for the 4 orbits in 500 s bins. Green regions denote high flux 
segments and red regions denote low flux segments. These high and low flux segments will be used for the lag analysis shown in Fig. [3] 
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Figure 3. Lag-frequency spectrum for the low (red) and high 
(green) flux segments shown in Fig. [2] The most negative lag for 
the high flux intervals is -660 ± 280 s at v = 1.8 X 10^"' Hz, and 
the most negative lag for the low flux is higher: —230 ± 31 s at 
u = 7 X 10"* Hz. Again, note that the most negative lag for the 
total observation (Fig. [1} occurs at p = 4.1 X 10"'' Hz. 



Figure 4. Lag-frequency spectrum of the 4th orbit alone, for the 
low (red) and high (green) flux segments shown in Fig. [2] The 
most negative lag for the high flux interval is —350 ± 130 s at 
u = 2.2 X lO"** Hz, and the most negative lag for the low flux is 
higher: -118 ± 42 s at = 6.5 X 10""' Hz. 



servation. In this section, we examine the lags from low and 
high flux intervals. Fig. [5] shows light curves from the 4 or- 
bits. The low flux sample is shown in red, and the high flux 
is shown in green. We chose light curve segments that were 
long in order to obtain information at low frequencies. The 
power spectral density confirms that up to 1.5 x 10""^ Hz, 
we are not in a Poisson noise dominated frequency regime 
for both the low and high flux samples. 



3.2.1 The lag -frequency spectrum 

Fig. [3] shows the lag vs. frequency for the high and low 
flux segments. The lag is computed in the same way as in 
Section 13.11 between the energy bands of 0.3 - 1 keV and 
1.2 - 5 keV. Focusing on the high flux points in green, we 
see that the soft lag occurs at lower frequencies than found 
from the total light curve sample, in Fig. [T] This means that 
the variability associated with the soft lag occurs on longer 
timescales during the flaring period. Also, the amplitude of 
the lag is greater. Now looking at the red low flux points, 
we see that the soft lag occurs at higher frequencies than 



the total, indicating the source variability is occurring on 
shorter timescales. 

As a check, we compute flux-dependent lags using orbit 
4 alone, which shows distinct low and high flux segments (as 
designated in green and red in the rightmost panel of Fig. [2)). 
The lag is computed for these two sets of continuous light 
curves, and shown in Fig. |4] We flnd that for continuous 
segments (i.e. for one realisation of some underlying pro- 
cess), the same flux-dependent behaviour of the lag is clear, 
just with lower signal-to-noise. The analysis that follows has 
been performed with low and high flux segments from the 
total 500 ks observation. 



3.2.2 The lag-energy spectrum 

We also examine how the flux-dependent lag evolves with en- 
ergy at the frequencies of the negative lag (i.e. [5.8 — 10.5] x 
10"* Hz for the low flux, and [1.4 - 2.8] x 10"" Hz for high 
fluxes). In this analysis, we measure the time lag of light 
curves in relatively narrow energy bins with respect to the 
ligh t curve of a broad er reference band from 0.3 - 0.8 keV 
(see lKara et al.|[2012l . for further details). We use the con- 
vention that a positive lag means that the light curve in that 
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Figure 5. Energy spectra for the high and low flux samples in panels (a) and (b) respectively. For clarity, only the power-law and 
reflection components are plotted. Panel (c) shows the powerlaw and reflection for the low and high flux spectra over plotted. 
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Figure 6. Lag-energy spectra for the low-flux (red) and high-flux 
(green) intervals, showing the energy dependence of the lag for 
the frequencies of the soft lag {u = [5.8 — 10.5] X 10"* Hz for the 
low flux, and u = [1.4 - 2.8] X 10"'' Hz for the high flux). 



energy bin lags behind the reference band, while a negative 
lag means that the Ught curve in that energy bin leads the 
reference. 

Fig. [6]shows the high-frequency lag-energy spectrum for 
the low and high flux segments. The low flux segments show 
a much clearer signal than the high flux, likely due to poor 
statistics and lower coherence in the high flux. The general 
trends of the two lag-energy spectra are similar, except that 
the amplitude of the lag between 1-4 keV is greater for the 
high flux. The low-flux lag-energy spectrum shows a clear 
peak at ~ 6.5 keV, the energy of the Fe Ka line. There 
does not appear to be a corresponding peak at high fluxes, 
but as the error bars are so large, it cannot be ruled out. 

3.2.3 The time-integrated energy spectrum 

We complement the lag analysis presented above with a brief 
investigation of the spectral properties in the low and high 
flux segment s . We a pplied the best fit spectral model from 
iFabian et all l|2012bl i to the low and high flux spectra, freez- 
ing all physical parameters that are not likely to change 
between state (i.e. inclination, galactic absorption, Fe abun- 
dance). The results are shown in panels (a) and (b) of Fig. [5] 



for the high and low flux segments, respectively. For clarity, 
only the power-law and total reflection components are plot- 
ted. The high flux spectrum is dominated by the powerlaw, 
while the low flux spectrum, shows a greater contribution 
from reflection. 

Panel (c) of Fig. [5] shows the power-law and reflection 
components of the low and high flux segments overplotted. 
We notice that the power-law flux increases by an order of 
magnitude between the low and high flux spectra, but the 
reflection component remains nearly unchanged, especially 
at the energy of the Fe K line. 

In lKara et al.l (|2012l ). we use Monte Carlo simulations to 
show that dilution will have an effect on the measurement of 
the lag. The lag is measured between a soft energy band and 
a hard band, where, to flrst order, we approximate the soft 
band as reflection and the hard band as powerlaw. However, 
we understand that there is 'contamination' in each band 
caused by the other varying component. As these compo- 
nents have different intrinsic time delays, the net measured 
lag will be some weighted average of all the components. 
In the ideal case, where there is no contamination, the in- 
trinsic lag will be the measured lag. In the other extreme, 
where there are equal parts of reflection and powerlaw in 
each band, the measured lag will be zero. Furthermore, con- 
tamination by a non-varying or an uncorrelated component 
will not effect the measurement of the lag, as the lag is mea- 
sured between coherent signals. 

It is clear from the spectra in Fig. [5] that we need to 
account for dilution in the measurement of the lag. Here we 
quantify the amount of dilution by measuring the reflection 
fraction in the hard and soft bands. For example, in the low 
flux spectrum, the soft band is composed of 80% reflection, 
but also the hard band is composed of 45% reflection. The 
net amount of reflection between the two bands is therefore 
35%, meaning that the measured lag is only 35% of the 
intrinsic lag. Similarly for the high flux spectrum, the soft 
band is composed of 40% reflection, while the hard band 
is contaminated by 10% reflection. Again the net amount of 
reflection is 30%, meaning that the measured lag in the high 
flux is only 30% of the intrinsic lag. We conclude then that 
the effect of dilution is nearly the same for both the high 
and low flux lags, and therefore we can be confident that 
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Figure 7. The lag-frequency spectrum for the total observation 
(same as Fig. [1} in red, o verplotted with th e lag-frequency spec- 
trum of 1H0707-495 from lKara et al.l 1I2OI2I) in blue. 



Figure 8. High-frequency lag-energy spectrum for the low flux 
segments in IRAS 13224-3809 (red), overplotted with the high- 
frequency lag-energy spectrum of 1H0707-495 {u = [0.98 — 2.98] X 
10-3 Hz) in blue. 



the intrinsic amplitude of the lag in the high flux sample is 
indeed greater than the lag in the low flux. 



4 DISCUSSIONS 

(i) As shown in iFabian erahl (|2012bh , IRAS 13224-3809 
exhibits frequency dependent lags where the hard band lags 
the soft band at low frequencies below ~ 2 x 10"'' Hz, 
while the soft band lags the hard band at higher frequencies 
(Fig. ID). 

(ii) In this work, we take a step further and study the 
flux-dependence of the observed lags. We find that the soft 
lag has a larger amplitude and occurs on longer timescales 
during high flux states than during low flux states (Fig. 

(iii) Accor d ing t o the lag inte rpretation presented by 
iFabian et al] l|2009l ). IZoghbi et akl pJlQ ) and others, soft 
lags are reverberation lags due to the light travel time 
between the primary source and the inner accretion disc. 
Therefore, the low flux state (smaller reverberation lags 
at shorter timescales) refers to a compact emitting source 
that is closer to the central black hole, while the high flux 
state (bigger lag at longer timescales) refers to an extended 
emitting region further from the black hole. Moreover, we 
point out that a more compact corona will inevitably ir- 
radiate smaller disc ra dii because of light bending (e.g. 
iMiniutti fc Fabianll2004l ). thus providing a qualitatively self- 
consistent explanation for both the higher frequency and 
smaller amplitude of the soft lags during low flux states. 

The spatial extent of the corona can be interpreted as an 
increase in the height of the corona above the disc, or some 
radial expansion of the corona, or some combination of the 
two. An increase in the amplitude of the lag is easily un- 
derstood for a vertically extended corona as the light travel 
between the corona and the disc increases. For a radially ex- 
tended source, one can imagine that the corona irradiates a 
larger portion of the disc, causing a longer average lag. Some 
complexities may arise with this interpretation because the 
emissivity decreases greatly with larger disc radius, however 
we cannot rule out a radially extended corona from the ob- 
servations presented here. Both the lag amplitude and the 
frequency change by a factor of ~ 3.5 between low and high 



fluxes, which could imply that the corona extends 3.5 times 
further away in the high state than in the low. 

The energy spectra of the low and high flux segments 
(Fig.[S]) agree with the interpretation of an extended corona. 
The low flux sample shows a stronger contribution from 
reflection. This is expected from a more compact corona, 
where light bending effects close to the central black hole 
cause a greater fraction of continuum emission to be directed 
towards the disc. Spectral modelling of 1H0707-495 during 
a low flux state in 2 011 also shows a more compact corona 
(|Fabian et al.ll2012al ). 

The spectral difference between low and high flux sam- 
ples appears to be caused by a change in powerlaw flux, and 
not from a change in reflection. Long timescale variations 
in continuum that are decoupled from reflection have also 
been observed in M CG-6-30-15 (|Vaughan fc FabianI [20o3 : 
IMiniutti et al.|[2007l ). and explained within the framework 
of strong light bending. We note that while the reflec- 
tion spectrum appears unchanged on these long timescales, 
the frequency- resolved covarian ce spectrum (as shown for 
1H0707-495 in iKara et"alll2012l ) does show that reflection 
varies with the continuum on short timescales. 



4.1 Comparison with 1H0707-495 

The lag analysis of IRAS 13224-3809 shows striking similari- 
ties to the results found for 1H0707-495 using 1.3 Ms of data 
(|Kara et al.ll20lj ). Fig. [7] shows the lag- frequency spectrum 
of IRAS 13224-3809 overplotted with 1H0707-495 (note dif- 
ferent y-axis scales). We flnd that IRAS 13224-3809 shows 
same general positive-to-negative trend as seen in 1H0707- 
495, just with a larger amplitude negative lag that occurs at 
lower frequencies. The negative lag in 1H0707-495 is mea- 
sured to be -31.9 ± 4.2 s at 1.33 x 10"^ Hz, while the neg- 
ative lag in IRAS 13224-3809 is a factor of ~ 3.33 times 
more (—92 ± 31 s) at a frequency that is ~ 3.25 times lower 
(4.1 X 10"* Hz). Both the lag and the frequency change 
roughly by a factor of 3.3. Assuming a simple mass scaling 
relationship, this implies that IRAS 13224-3809 is 3.3 times 
more massiv e that 1H0707-495. Ac cording to the mass scal- 
ing trends in lDe Marco et al.l (|2012l ). the lags in IRAS 13224- 
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3809 are consistent with a black hole mass of 10 Mq . In this 
work, however, we use IRAS 13224-3809 to show that there 
is some flux-dependence to the reverberation lag. Therefore, 
we should take some caution in predicting the mass of the 
central black hole from the amplitude and frequency of the 
lag. While observations used to measure reverberation lags 
are relatively long, it is possible to sample some intrinsically 
high or low flux state. This would introduce systematics to 
the simple mass scaling relationship. 

We compare the low-flux lag-energy spectrum of 
IRAS 13224-3809 with the lag-energy spectrum of 1H0707- 
495 (Fig.|5}. The shapes of the spectra for these two different 
sources are strikingly similar, but the amplitude of the lags 
is greater for the low flux segments in IRAS 13224-3809. 
The peak at ~ 6 keV and the dip at 3 - 4 keV are clearly 
present in both data sets. Being able to isolate the low- 
flux segments of IRAS 13224-3809 has allowed us to clearly 
uncover the reverberation signal, as was only possible with 
1.3 Ms of data with 1H0707-495. 1H0707-495 does not have 
such distinct low-flux segments as IRAS 13224-3809, and 
so a similar flux-resolved analysis was not possible for this 
source. 

The high-frequency lag-energy spectra of IRAS 13224- 
3809 and 1H0707-495 (Figs. |6] and [H]! are consistent with a 
reflection model, where the negative lags are associated with 
the light travel time between the continuum-emitting corona 
and the accretion disc. The lag-energy spectrum shows the 
1-4 keV band (powerlaw dominated) leads the 0.3 - 0.8 
keV reference band (where reflection is appreciable) . Similar 
to 1H0707-495, we see a distinct peak at ~ 6 keV, the rest 
frame energy of the Fe K line, and a dip at 3 - 4 keV, which 
has been interpreted as a signature from the red wing of the 
Fe K line, originating from the innermost radius. This inter- 
pretation is discussed further in lKara et al.l (|2012l ). where we 
model the high-frequency lag-energy spectrum of 1II0707- 
495. 



4.2 Conclusions 

Looking at the flux-resolved reverberation lags reveals new 
observational constraints that can help us distinguish be- 
tween the physical nature of the low and high flux states. 
We learn that flares are likely connected to the expansion 
of the corona out vertically and possibly radially, too, and 
that when the source is in quiescence, the corona remains 
compact and close to the central region. 



ACKNOWLEDGEMENTS 

This work is based on observations obtained with XMM- 
Newton, an ESA science mission with instruments and 
contributions directly funded by ESA Member States and 
NASA. EK thanks the Gates Cambridge Scholarship. ACF 
thanks the Royal Society. 



Boiler T., Tanaka Y., Fabian A., Brandt W. N., Gallo L., 
Anabuki N., Haba Y., Vaughan S., 2003, MNRAS, 343, 
L89 

Chartas G., Kochanek C. S., Dai X., Poindexter S., 

Garmire G., 2009, ApJ, 693, 174 
Chartas G., Kochanek C. S., Dai X., Moore D., Mosquera 

A. M., Blackburne J. A., 2012, ApJ, 757, 137 
Dai X., Kochanek C. S., Chartas G., Kozlowski S., Morgan 

C. W., Garmire G., Agol E., 2010, ApJ, 709, 278 
De Marco B., Ponti G., Uttley P., Cappi M., Dadina M., 

Fabian A. C, Miniutti G., 2011, MNRAS, 417, L98 
De Marco B., Ponti G., Cappi M., Dadina M., Uttley 

P., Cackett E. M., Fabian A. C, Miniutti G., 2012, 

'arXiv:1201.0196 

Emmanoulopoulos D., McHardy I. M., Papadakis I. E., 

2011, MNRAS, 416, L94 

Fabian A. C, et al., 2009, Nature, 459, 540 

Fabian A. C, et al., 2012, MNRAS, 419, 116 

Fabian A. C, et al., 2012, arXiv:1208.5898vl 

Gallo L. C, Boiler T., Tanaka Y., Fabian A. C, Brandt 

W. N., Welsh W. F., Anabuki N., Haba Y., 2004, MNRAS, 

347, 269 

Guilbert P. W., Rees M. J., 1988, MNRAS, 233, 475 
Haardt F., Maraschi L., 1991, ApJ, 380, L51 
Haardt F., Maraschi L., 1993, ApJ, 413, 507 
Jansen F., et al., 2001, A&A, 365, LI 

Kara E., Fabian A. C, Cackett E. M., Steiner J. F., Uttley 
P., Wilkins D. R., Zoghbi A., 2012, MNRAS, 264 

Kotov O., Churazov E., Gilfanov M., 2001, MNRAS, 327, 
799 

Lightman A. P., White T. R., 1988, ApJ, 335, 57 
McHardy I., 1988, MmSAI, 59, 239 

McHardy I. M., Papadakis I. E., Uttley P., Page M. J., 

Mason K. O., 2004, MNRAS, 348, 783 
Miyamoto S., Kitamoto S., Mitsuda K., Dotani T., 1988, 

Natur, 336, 450 
Miniutti G., Fabian A. C, 2004, MNRAS, 349, 1435 
Miniutti G., et al., 2007, PASJ, 59, 315 
Nowak M. A., Vaughan B. A., 1996, MNRAS, 280, 227 
Nowak M. A., Vaughan B. A., Wilms J., Dove J. B., Begel- 

man M. C, 1999, ApJ, 510, 874 
Papadakis I. E., Nandra K., Kazanas D., 2001, ApJ, 554, 

L133 

Ponti G., et al., 2010, MNRAS, 406, 2591 

Ross R. R., Fabian A. C, 2005, MNRAS, 358, 211 

Shakura N. I., Sunyaev R. A., 1973, A&A, 24, 337 

Striider L., et al., 2001, A&A, 365, L18 

TerreU J., 1967, ApJ, 147, 827 

Uttley P., Wilkinson T., Cassatella P., Wilms J., 
Pottschmidt K., Hanke M., Bock M., 2011, MNRAS, 414, 
L60 

Vaughan S., Fabian A. C, 2004, MNRAS, 348, 1415 
Zoghbi A., Fabian A. C, Uttley P., Miniutti G., Gallo L. C, 

Reynolds C. S., Miller J. M., Ponti G., 2010, MNRAS, 401, 

2419 

Zoghbi A., Fabian A. C, Reynolds C. S., Cackett E. M., 

2012, MNRAS, 422, 129 



REFERENCES 

Arevalo P., Uttley P., 2006, MNRAS, 367, 801 

Boiler T., Brandt W. N., Fink H., 1996, A&A, 305, 53 



© 2012 RAS, MNRAS 000, [He] 



